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^_^ Abstract. By modelling quantum systems as emerging from a (classical) sub-quantum 

r^ thermodynamics, the quantum mechanical "decay of the wave packet" is shown to simply 

^ I result from sub-quantum diffusion with a specific diffusion coefficient varying in time due to 

J_^ a particle's changing thermal environment. It is thereby proven that free quantum motion 

JT^ strictly equals ballistic diffusion. The exact quantum mechanical trajectory distributions and 

Kxv the velocity field of the Gaussian wave packet are thus derived solely from classical physics. 

• Moreover, also quantum motion in a linear (e.g., gravitational) potential is shown to equal said 

fj ballistic diffusion. Quantitative statements on the trajectories' characteristic behaviours are 

• i-H obtained which provide a detailed "micro-causal" explanation in full accordance with momentum 

K^_^ conservation. 

►^ 1. Quantum mechanical dispersion of a free Gaussian ^vave packet exactly 

QQ modelled by sub-quantum ballistic diffusion 

•/"J Considering quantum theory as emergent [1], we propose that it may result from a deeper, 

more exact theory on a sub-quantum level. In our approach, one assumes that the latter can 
be described with the aid of non-equilibrium thermodynamics. We ask ourselves how quantum 
theory would have evolved, had the "tool" of modern non-equilibrium thermodynamics existed, 

^ say, a century ago. As one of us (G. G.) has recently shown, one can derive the exact Schrodinger 

equation with said tool, where the relation between energy E and frequency oj is used as the 
only empirical input, E = hu ^^, with the additional option that even the appearance of 
Planck's constant h may have its origin in classical physics [1]. For an extensive review of refs. 
[2] and [3], for connections to similar work, and, in particular, to Fisher information techniques, 
C^ see [I]. As to approaches in a similar spirit, see, for example, O El [71 El EJ HOl IH [E]- 

In the thermodynamic approach to quantum behaviour [21 El H] , a particle of energy E = Hw 
is characterized by an oscillator of angular frequency cj, which itself is a dissipative system 
maintained in a non-equilibrium steady-state by a permanent troughput of energy, or heat flow, 
respectively. The latter is a form of kinetic energy different from the "ordinary" kinetic energy 
of the particle, as it represents an additional, stochastic contribution to it, i.e., from the presence 
of zero point fluctuations. Therefore, one deals with a context-dependent total energy of the 
whole system (i.e., the particle as the "system of interest" in a narrower sense and the heat flow 



constituting the particle's noisy thermal embedding), which is assmned as 

i?tot = ^ + ^, (1) 

where 6p := mu is said additional, environment-related fluctuating momentum of the particle 
of mass m. 

We let ourselves be inspired by the experiments with the classical "walkers" of Yves Couder's 
group (see, for example, [13l[lll[15l[16]). One has to do with a rapidly "bouncing" (oscillating) 
object, which itself is guided by an environment that contributes some fluctuating momentum 
to the walker's propagation: a bouncing droplet phase-locked with an oscillating fluid surface. 
However, the walker is also the cause of waves emanating from the particle "hitting" the 
surface, and it is the detailed structure of the emerging wave configurations, in turn, that 
stochastically influences the walker's path. Thus, if we imagine the bouncing of a walker in its 
fluid environment, the latter will become "excited" or "heated up" wherever the momentum 
fluctuations direct the particle to. After some time span, a whole area of the particle's 
environment will be coherently "heated up" in this way. We shall here try to implement a 
similar "particle-wave coupling" with respect to a quantum system. 

To start with a first major aspect of our modeling, we consider a particle of initial velocity 
V in the midst of "well ordered" diffusion waves, which themselves emerge out of the noisy, 
Brownian-type diffusions of myriads of single sub-quantum particles. Being swept along with 
a diffusion wave, with initial {t = 0) location x(0) and diffusion velocity u, a thus described 
"quantum" particle's distance to some heat accumulation's center xq at time t will be 

x{t) = x{0) + ut , (2) 

such that one obtains the r.m.s. of ([2| as 



a;2(i) = a;2(o) + 2 x(0) • u{x, t) t + u'^{x, t) r . (3) 

Now we introduce an essential point: we assume, as an emerging result out of the statistics 
of a vast number of diffusion processes, the complete statistical independence at any time t of 
the velocities u and the positions x, and thus also of u and v. 

xu = and vu = . (4) 

This is justified considering the statistics of huge numbers, millions of millions of diffusive sub- 
quantum Brownian motions, which are supposed to bring forth the emergence of said larger- 
scale collective phenomenon, i.e., the diffusion wave fields as solutions to the heat equation [3j. 
Therefore, with the average orthogonality of classical (convective) momentum mv on one hand, 
and its associated diffusive momentum mu on the other, one gets rid of the term linear in t in 
Eq. ([3]) , and thus of irreversibility, and one obtains 



a;2(i)=a;2(0)+n2t^ (5) 

Eq. ([5]) is the result obtained for the "pure" emergent diffusive motion as given by ([2]). 

Now, let us assume we have a source of identical particles, which are prepared in such a way 
that each one ideally has an initial (classical) velocity v. Even if we let them emerge one at a 
time only, say, from an aperture with unsharp edges (thus avoiding diffraction effects to good 
approximation), the probability density P will be a Gaussian one. This comes along with a heat 
distribution generated by the oscillating ("bouncing") particle(s), with a maximum at the center 



of the aperture xq = vt. As the classical diffusion equation, in one dimension for simplicity, is 
given by 

d 52 

where D is the "diffusion constant", we have the corresponding solution 



-Pix,t) = D—Pix,t), (6) 



1 (^-^0) 

P(x,t) = -^^e ^^^^ (7) 

V2Tra 

with the usual variance a^ = x^(t). Note that from Eq. ([l]) one has for the averages over particle 
positions and fluctuations (as represented via the probability density P) 



{6p) 



2 



Etot = huj + — — = const., (8) 

2m 

with the mean values (generally defined in n— dimensional configuration space) 

JSpf := I P{6pf d"x . (9) 

Equation (Is]) is a statement of total average energy conservation. This means that a variation 
in 5p implies a varying "particle energy" fujj, and vice versa, such that each of the summands 
on the right hand side for itself is not conserved. In fact, as shall be detailed below, there will 
generally be an exchange of momentum between the two terms providing a net balance 

m5v = mSu , (10) 

where 6v describes a change in the "convective" velocity v paralleled by the "diffusive" 
momentum fluctuation 6{5p) := m6u in the thermal environment. 

As elaborated in references [2] and [3], once Eq. (II]) is assumed, considerations based on 
Boltzmann's relation between action and angular frequency of an oscillator provide, without 
any further reference to quantum theory, that 

5p = mu := hka = VlnP . (11) 

Further, as can easily be checked by integration, it holds that 



(VlnP)2 = -V^lnP. (12) 



As in Eq. ([8]) only the kinetic energy varies, one has 5E]^\^{t) = (5£'km(0) = const. Thus, with 
the Gaussian M, this provides an expression for the averaged fluctuating kinetic energy, or heat, 
of a particle and its surroundings, 



m , ^ ,„ m—F: m- 



<5Skin(i) = -wiSvf + -u^ = -[5vf + 
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5^kin(0)= 0+1^.2 
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Moreover, with the diffusion constant 

D:=— (14) 

2m ^ ' 



Eq. ( 13 ) provides an expression for the initial vefocity fluctuation 



UQ = — . (15) 

However, in a next step we now take into account the small momentum fluctuations mJn, 
providing an altered convective velocity v — )• v ^ bvit)^ and thus an additional displacement 
bx = \6u\ t = \6v\ t, i.e., as soon as t > 0. Therefore, returning to the simple rule of Eq. ([2]), 
one now must decompose u{t) into its initial value uq and a fluctuating contribution 5u{t), 
respectively. Unless some thermal equilibrium were reached, the latter is typically given off 
from the "heated" thermal bath to the particle of velocity v, 

u{t) = uo- 6u{t) . (16) 

As is detailed in ref. [IJ], this leads ultimately to the substitution of u by uq in Eq. ([5]), i.e.. 



x^{t) = x^{0) + u'ot' . (17) 



Inserting (15) into (17) for the particular case that x'^{t) = a"^ (i.e., x'^{0) = (Jq) provides for the 



time evolution of the wave packet's variance 



'^'o 



a^ = aUl + ^\ . (18) 



The quadratic time-dependence of the variance a^ is remarkable insofar as in ordinary diffusion 
processes the scenario is different. There, with the Gaussian distribution being a solution of the 
heat equation, for purely Brownian motion the variance grows only linearly with time, i.e., as 
described by the familiar relation 



x^{t) = x^{0) + 2Dt . (19) 

However, as we have seen, the momentum exchange between the particle and its environment 
is characterized by both a changing velocity and by a changing thermal environment of the 



particle, i.e., also by a changing diffusivity. Therefore, Eq. (19) must be modified to allow for 
a time-dependent diffusivity. In other words, we shall have to deal with the field of anomalous 
diffusion. (For a short introduction, see, for example, ref. pLSj.) This means that instead of 
the diffusion constant D in Eq. (|6]), we now introduce a time-dependent diffusion coefficient 
D{t) = kt°^, where A; is a constant factor. The exponent a was derived [TU] to be a = 1. 
Therefore, k = -D^/cTq = Uq, and the time-dependent diffusion coefficient becomes 

D(t)^4t=^t = J^^t. (20) 

Note that with the exponent of t being a = 1, or the t^-dependence of a^ in (18), respectively. 



one deals with the special case of anomalous diffusion usually named ballistic diffusion. We 
shall review some general properties of ballistic diffusion at the end of the paper. At this point, 
however, it is useful to recall that throughout our modelling of sub-quantum processes, we 
deal with various processes at different time scales. On the shortest scales, we have assumed 
Brownian- type motions (not detailed here), which, on the next higher level of (spatial and) 
temporal scales lead collectively to the emergence of a regular diffusion wave. The latter is 
characterized by a velocity u according to ^, and it is orthogonal on average to the particle's 
velocity v, thus providing the r.m.s. displacement ([S]) depending on u{t). As a next step, we 



have introduced the noisy thermal bath of the particle's environment, i.e., essentially the effect of 
other diffusion wave configurations, which disturbs the relation ([5| by introducing a fluctuating 



term 5u. The net effect of the latter, however, is the r.m.s. displacement (17) with a dependence 
solely on the initial diffusive velocity uq. This manifests itself also in the expression for D(t) of 
the ultimately emerging ballistic diffusion, which is also dependent only on uq. However, even 
on this level of ballistic diffusion one can recover the signature of Brownian motion. In fact, if 
one considers the time-average of D{t) for large enough times t ^ 1/w, i.e., 

{D{t)):=]I^D{t')dt'='§t=^, (21) 

one immediately obtains the linear-in-time Brownian relation 

x2(f) = x2(0) + 2(D(t)) t and a^ = a^ + 2{D{t))t , (22) 



which is, however, also in accordance with the t^-dependence of Eq. (17). 
Note that the diffusivity's rate of change is a constant, 

dD(t) D^ 2 

— - — = —2" = Uq = const., (23) 

such that it is determined only by the initial r.m.s. distribution ctq, providing also a re-formulation 



of Eq. (17) as 



dD{t) 



2;2(t)=2;2(0) + ^^^^^t2_ (24) 

dt 



With the square root of (18), 



a = ao\ll + ^ , (25) 



^0 



we note that cr/ao is a spreading ratio for the wave packet independent of x. This functional 
relationship is thus not only valid for the particular point x{t) = cr{t), but for all x of the 



Gaussian. Therefore, one can generalize (25) for all x, i.e.. 



x(t) = x(0) — , where — = J 1 + ^^ . (26) 

Now we remind ourselves that we deal with a particle of velocity v = p/m immersed in a 
wave-like thermal bath that permanently provides some momentum fluctuations Sp. The latter 



are reflected in Eq. (25) via the r.m.s. deviation a{t) from the usual classical path. In other 



words, one has to do with a wave packet with an overall uniform motion given by v, where the 



position xo = vt moves like a free classical particle. As the packet spreads according to Eq. (25), 
x{t) = a{t) describes the motion of a point of this packet that was initially at x(0) = ctq. 
Depending on whether initially x(0) > ctq or x(0) < (Tq, then, respectively, said spreading 
happens faster or slower than that for x(0) = o"o. In our picture, this is easy to understand. For 
a particle exactly at the center xq of the packet [x(0) = 0], the momentum contributions from 
the "heated up" environment on average cancel each other for symmetry reasons. 

However, the further off a particle is from that center, the stronger this symmetry will be 
broken, i.e., leading to a position-dependent net acceleration or deceleration, respectively, or, in 
effect, to the "decay of the wave packet". Moreover, also the appearance of the time-dependent 
diffusivity D{t) is straightforward in our model. Essentially, the "decay of the wave packet" 
simply results from sub-quantum diffusion with a diffusivity varying in time due to the particle's 



changing thermal environment: as the heat initially concentrated in a narrow spatial domain 
gets gradually dispersed, so must the diffusivity of the medium change accordingly. 

Finally, one obtains with Eqs. (26) and (14) for the "smoothed out" trajectories (i.e., those 



averaged over a very large number of Brownian motions) 



XtotW = vt + x{t) = vt + x{0)— = vt + x{0)Jl + —^ 



4 • 



(To y 4?T1,^(Tq 

Also, one can now calculate the average total velocity of a Gaussian wave packet, 

.,x dj;tot(t) f.s . dx{t) 
v^t) = -^^ = v{t) + -^ , 



(27) 



providing [ITj 



Wtot(t) = v{t) + [Xtot(i) - Vt] 



h^ t 



Am? ct^iTq 



(28) 



(29) 



2. Simulation of quantum phenomena w^ith purely classical diffusion using coupled 
map lattices 



It is straightforward to simulate the diffusion process of Eq. ( 22 ) in a simple computer model. 



Using coupled map lattices (CML), one approximates the heat equation as usual by 
P[i, k + l]= P[i, k] + ^^J^ {P[i + 1, fc] - 2P[i, k] + P[i - 1, k]} , 



(30) 



with space and time grid indices i and k, respectively. For our anomalous ("ballistic") diffusion 



one simply inserts (20) into (30). 

The result is depicted in Fig. [II where the (macroscopic, classical) velocity is chosen as 
V = 0. (For examples with v ^ and different co, see [IZ].) Moreover, nine exemplary 
averaged Bohmian trajectories are shown in Fig. [T| and it must be stressed that the Figure 
shows the emerging behaviour of the Gaussian packet following solely from the CML simulation 



of Eq. (30). In addition, the emerging trajectories from the simulation (dark lines) are shown 



together with the calculated ones from Eq. (27) (bright, mostly superposed by the dark lines). 



providing exactly the same trajectories (i.e., up to resolution limits due to discretization). 

Note that the trajectories are not the "real" ones, but only represent the averaged behaviour 
of a statistical ensemble. The results are in full concordance with quantum theory, and in 
particular with Bohmian trajectories. (For a comparison with the latter, see, for example, [20], 
or the Figures for the Gaussian wave packet example in [21], which are in excellent agreement 
with our Fig. [!}) This is so despite the fact that no quantum mechanics has been used yet, 
i.e., neither a quantum mechanical wave function, or the Schrodinger equation, respectively, 
nor a guiding wave equation, nor a quantum potential. Moreover, we want to stress that our 
model offers possible insights into the sub-quantum domain which must escape (Bohmian or 
orthodox) quantum theory because the latter simply does not employ the "language" necessary 
to express them. Note, for example, that the existence of the hyperbolic trajectories depicted 
in Fig. q1 which are given by the formula for the scale invariant wave packet spread (18), has 
a simple physical explanation in terms of sub-quantum processes. As the inflection points of 
the hyperbolas are, according to (18), characterized by the relation D'^t^/a^ = 1, i.e., by the 



length scales uit = Uq, the trajectories' evolution is easily understood: as long as the main bulk 



0' 



of the heat "stored" in the initial Gaussian spreads well "inside" the distribution, Ugt < a, 
the average particle velocity v is not affected much. However, if said main bulk approximately 
reaches the distance ctq, or spreads to regions Ugt^ > ctq, respectively, the particles will "feel" the 





Figure 1. Classical CML simulation of 
the dispersion of a free Gaussian wave 
packet. 



Figure 2. Dispersion of a Gaussian 
wave packet in a linear potential (e.g., a 
gravitational field). 



full heat and get propagated into new directions. For t 
rate of the whole Gaussian packet: 



CO, then, uq becomes the spreading 



da 
~dt 



hH 



i— >oo 
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Uq . 



(31) 



In other words, the "spreading" already begins at i = 0, but becomes "visible" in terms of 
deflected trajectories only when t = ctq/uq. 

So far we have shown that free one-particle quantum motion is identical to sub-quantum 
ballistic diffusion. This is the basis of a research program that would eventually cover more and 
more complex situations beyond the case of free motion. As a first simple example, we extend 
the present scheme to include a linear potential. That is, we place the initial Gaussian packet ^ 
in a uniform potential V = K ■ x, which may be an electric or a gravitational field, for example. 
For illustration, but without loss of generality, we substitute in the following K by mg, i.e., we 
deal with a Gaussian packet freely falling due to the potential V = mg ■ x. 

After a short calculation p^, one obtains the trajectories of particles in a gravitational field. 



Xtot{t) = vt 



9-2 . ^ / ^^i^ 



(32) 



and the particle acceleration 

a;tot = -g + x(0)-2 

^0 



1 + 



^ 

-i 



-3/2 

(33) 



In Fig. [2j exemplary trajectories of a Gaussian in a gravitational field are shown as obtained 
by the CML simulation of ballistic diffusion, modified by the substitution v ^ v — gt. The 



trajectories (dark lines) superpose those (bright lines) derived directly from Eq. (32), and again 
exhibit excellent agreement. 

Note that some trajectories of the dispersing Gaussian even overcome "gravity" for a well- 
defined period of time, as can also be seen in Fig. [2} In fact, our sub-quantum model provides 
a detailed explanation of why, and within which time limits, this "anti-gravity" effect becomes 



possible. A look at the last term of Eq. (33) provides the answer. Similarly to the discussion of 
the hyperbolas' inflection points in the free case, one deals also here with an extremum at the 
scale UqI'^ = Uq. However, this time the corresponding expression (in rectangular brackets) is 
antagonistic to g. In other words, said scale describes the maximum of the "anti-gravity" effect, 
because it is there where the heat of the main bulk of the packet is consumed, which has via 
the kinetic energy counter-acted the effect of gravity. For larger times, then, the remaining heat 
gets gradually less, and therefore gravitational acceleration begins to dominate the trajectories' 
curvature. 

Finally, instead of following just one Gaussian, we extend our simulation scheme to include 
two possible paths of a particle which eventually cross each other. For this, we use two Gaussians 
approaching each other. In the space-time diagram of Fig. 3, they are characterized by the same 
half-widths a, whereas in Fig. |4]we choose Gaussians whose half-widths' ratio is 2:1. 

Note that in both Figs . p^ and |4] one can observe a basic characteristic of the (averaged) particle 
trajectories, which, just because of the averaging, are identical with the Bohmian trajectories. 
To fully appreciate this characteristic, we remind the reader of the severe criticism of Bohmian 
trajectories as put forward by Scully and others (see |22], and references therein.) The critics 
claimed that Bohmian trajectories would have to be described as "surreal" ones because of their 
apparent violation of momentum conservation. In fact, due to the "no crossing rule" for Bohmian 
trajectories in Young's double slit experiment, for example, the particles coming from, say, the 
right slit (and expected at the left part of the screen if "classical" momentum conservation 
should hold) actually arrive at the right part of the screen (and vice versa for the other slit). 
In Bohmian theory, this "no crossing rule" is due to the action of the non-classical quantum 
potential, such that, once the existence of a quantum potential is accepted, no contradiction 
arises and the trajectories may be considered "real" instead of "surreal" . 

Here we can note that in our sub-quantum approach an explanation of the "no crossing rule" 
is even more straightforward and actually a consequence of a detailed microscopic momentum 
conservation. As can be seen particularly clearly in Fig. ^ the (Bohmian) trajectories are 
repelled from the central symmetry line. However, in our case this is only implicitly due to a 
"quantum potential" , but actually due to the identification of the latter with a kinetic (rather 
than a potential) energy: it is the "heat of the compressed vacuum" that accumulates along 
said symmetry line (i.e., as reservoir of "outward" oriented kinetic energy) and therefore repels 
the trajectories. In Fig. 4, moreover, one sees the resulting trajectories for the case that the 
symmetry of Fig. 3 is broken: by using different half-widths for the two paths, one obtains the 
skew-symmetric distribution of particle intensities and trajectories, respectively. Both Figs. [3] 
and H^ are in full concordance with Bohmian trajectories (cf., for example, [23] for comparison). 
However, as mentioned, in our case also a micro-causal explanation is provided, which brings 
the whole process into perfect agreement with momentum conservation on a more "microscopic" 
level. 





Figure 3. Trajectories of two Gaussians 

with the same half widths. 



Figure 4. Trajectories of two Gaussians 
with half width ratio of 2:1. 



To sum up, we have shown for the cases of free motion and motion in linear potentials, 
respectively, that the time evolution of a one-particle quantum system in the noisy heat bath 
of the surrounding "vacuum" exactly equals that of (classical) ballistic sub-quantum diffusion. 
Note that there are some well-known general characteristics of ballistic diffusion [24J, and the 
results presented in this paper agree perfectly with them. For one thing, ballistic diffusion is the 
only type of diffusion that exhibits reversibility, and because of this it violates ergodicity (i.e., 
as in our cases). Also, if the ballistic system is not in equilibrium initially, it will never reach 
equilibrium (which is true here as well) . Finally, the result of any measurement depends on the 
initial conditions. This can be clearly seen also from our results for the time evolution of the 
Gaussians and the corresponding averaged trajectories, which all depend on the initial values of 
no, or (To, respectively. 
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